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THEORY OF AN AIRPLANE ENCOUNTERING GUSTS, ni. 
By Edwin Bidwecl Wilson. 



The following is an immediate continuation of my two previous papers on ths same sub- 
ject which appeared in these Reports, First Tear, pages 52-75, and Iliird Year, pages 406-431, 
the latter being a reprint of my contribution to the Ftoceedings American Philosophical Society 
(Philadelphia), volume 58, pages 212-248. In my second paper I pointed out: (1) That the 
study of a gust of the type Je^ sin 'pt, tuned both in damping and in period to the natxiral 
motion of the machine, might be important; (2) that the solution for this case by Brodetsliy 
treated only the particular integral without taking account of the fact that the constants of 
integration in the complementary function might ba such as largely to upset any conclusion 
that such a gust produces violent motions; (3) that a new melihod for solving linear equations 
had been developed by Bromwich, which was suited to determine the motion for any particular 
gust, when the machine started from equilibrium, without the trouble of determining the 
constants of integration in the complementary function. 

I shall now apply Bromwich's method to the calculation of the eJffect of a head-on gust of 
the form u.^ = Jr^ sin ft where n = .06S4, p=.187, as in the case of the slow oscillation for 
the Curtiss JN'2. Inasmuch as the method depends on tiie use of the theory of functions of a 
complex variable, which is a mathematical subject of prime importance to any aeronautical 
engineer who would have that knowledge of fluid motion which is regarded in high quarters as 
essential, and further inasmuch as neither the theory of functions nor Bromwich's special 
method is likely to become as familiar as they should be to engineers or physicists without 
detailed directions for and examples of the application of such ways of calculating, I may be 
pardoned for the somewhat lengthy presentation of matters elementary for tiie pure matiie- 
matician. 

Suppose it be required to solve the equation 

^+ 2u^ + (j^+n*)x = Fe-* cos (7rf+«) (1) 

with the supposition that the damped harmonic force on tiie right is applied at and from the 
time t = 0, and that at j=0 the system is at rest in its position of equilibrium, i. e., 0 and 
dx/dt = 0. 

Now trigonometric terms are generally treated by their exponential equivalents, through 
the formulas 

003 y=' — g — » Bmy=- — gj — > 

e*»=cos y+i sin y, where 

The solution of (1) may be obtained by solving the equation obtained by replacing 
cos (nt+a) by its value as the sum of two e:q>onential expressions. A method generally better 
is to replace the equation by 

§-t-2u§+ («» +jj»)a!- ^fe-*e*e^> (2) 
and take the real part of the solution since the real part of e*" is cos y. 
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Bromwich states that the solution of (2) subject to the stated conditions is 
where 

^°"(X»+2i;X+u*+»*)(X+u-in) 

■'(N+u-in)*(X+u+tn)' * 

and where the integral must be taken around a curve in the complex X-plane sufiGiciently large 
to include all the points, X = — u— in, and X = — u + in (double), where f becomes infinite. He 
further points out that the value obtained for x is none other than the sum of the residues of 
by virtue of Cauchy's Uieorem. 

Now a residue may be given a simple definition. Suppose 

where ^Oc) does not vanish when x=-a, i. e., no factor x — a may be canceled out, the fraction is 
in its lowest-terms. The function <^(z) may be expanded by Taylor's theorem about i = a as; 



(o) (x - a) -k-M" (a) Oc - a)* + 

and 

where 4<=<^iW(<i)/i! with 4>^^{a) denoting the ith derivative of 4>ix) taken for the value t=a. 
The coefficient which occurs over the factor z-a is called the "residue of /(a) at x—a." 

In case n = l, that is, M f{x) = the residue of/ at a is i^>(a). 

K^ — Q,) 

In case n=2, that is, if fix) °° (x-^^)* ' ^^'^^^ of / at a is ^' (a). 
As applied to the case in hand where 

fi\) ={«^*= (X + „ -in)\-K+v+iny 
the function for the consideration of the value \=—v—in may be written 



/(X)-- 



X + v+tn 



and the residue of / at X — —u—i» is obtained by substituting X = — u— in in the numerator; 
hence residue of /at i'W~"~*""{—2in)~*; and for the consideration of the value 

X— — u+m, 

' 

and the residue of/at X» — u+tn is the value for this value of X of the derivative: 

d Feh'^ ^ F^[{\+v+in)tf^-(^^ 
dk X+u+in"" (X-f-u+in)* 

which value is 

residue of /(X) at -w+t» = Fef-e'-»*"»«(2in<-l)(2ift)-*, 
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Hence the solution for x in (2) is the sum of the residues, or 

ar— ^«-' «^*-*"*+e*"*(2i7rf-l)]- 

And the solution for j in (1 ) is the real part of the above obt-xinea by substituting for «'•*, 
e*^ their expressions in terms of trigonometric functions. This gives: 



^^^nf sin (rat sin a sin nfj. 



the result stated by Bromwich, without giving the steps in detail. 
If it had been required to solve the equation: 

^+2i;^+ (u»+ ft*) »=- Fer'*taiL (ra+ a), (la) 

the only change would have been to take the imaginary part, instead of the real part, and 
throw out the factor i =••/— 1- Thus, 

-F 

*" 4^ e*- [2t7rf «<"«- ie*'*-€r*^] 
— F 

= (cos w+t sin ») [2tnt cos nt— 2n* sin n#— 2* sin nf] 

■= «^ (— 2iit cos M sin nt—27U sin w cos ni+2 sin u sin nt) + 

i (2ni cos a cos nt— 2«f sin u sin n< — 2 sin a sin n<) 
and the solution of (3) is 

F 

a;=25^ «~* [— fii cos (.nt+a) +cos a sin n£[. 
For the head gust the equations to be solved are 

- M^u- M„w + Q?jy* - Mfi) 9 - Jfj^ff-* sin pt. 

These are replaced by equations with Jer^+P* instead of Je-^ i* and the equationr 
for n, t become 

Q.-X^ i-X„r,-(X^+g) i-=X„//(X+ft-2») 
-Z„f +(X-Z„)i}- (Z5+ D) >f-Z,J/(\+ft-2n) 

- if„f- + (X;',x»- Ji;^) r= if;,j/(x+»-2n). 

Next If, f are obtained by solution and multiplied by The resulfa are as follows: 

„ - (.128X* - 1 .16X' - 3.386X- .917) Je" 

'°'(X*-8.49X«-24.6X»-3.385X-.917) (X+n-pi) ' 

„ -(X'.557X- 2.458) Je'^ 

°" (X* - 8.49X» - 24.5X* - 3.386X- .917) {X+ n- pi) ' 

_J -.02851Xjg» 

" TX*- 8.49X*- 24.5X»-3.385X- .917) (X+»- pi) ' 

The solutions for if, ff, respectively, are the sums of the residues of these three expressions 
The denominator factors (since »=.0654, p=-.187) into (X+ .0654+. lS7i)(X+. 0664- 
.lS7t)» (X-4.18+2.43i)(X+4.18— 2.43i). It is.necessary to calculate the residues for each o, 
the following values of X: - 4.18 ±2.43 i, - .0654 - .187i, - .0654 + .187i. The first three corre- 
spond to single factors of the denominator and are obtained merely by discarding that factor 
and substituting the value of X in the rest of the expression; the fourth requires that the double 
factor be discarded and that the value of X be substituted in the derivative of what remains. 
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The whole of the calculation need not be made. The interest lies almost entirely in the up 
and down motion which is given by 

and for which w and 8 alone need to be determined. One of the advantages of the Bromwich 
method lies precisely in this ability to calculate just those elements needed. 
We shall begin with 9 and figure the residues of: 

-.02851 XJe^ _ 

(X + .0664 + .187i) (X + .0664 - .187i)» (X + 4.18 ±2.43i) ' 

For that at -4.18+2.43 i calculate 

» -02851 (-4.18 + 2.43i) g-^-"* (cos 2.43t+i sin 2.43f) 
^«"" (-4.115+2.617i) (-4.116.+2.243i)» (4.86i) ' 

discarding for the moment the factor —J common to all residues. The rules for dealing 
with imiaginaries by trigonometry are helpful. We need the magnitudes and angles of the 
quantities: 

mag (-4.18+ 2.43i) -4.834, log mag =0.6843, ang=149*.83, 
mag (-4.115 +2.617i) =-4.875, log mag= 0.6879, ang=« 147''.54, 
mag (-4.115 +2 J!43i)- 4.792, log mag =0.6805, ang=151''.40. 

The angle of the coe£5cient in is -30''.51 ; the log mag is 6.4038. Hence 
5, = (.0002183 -.0001283 i) e^'Kcos 2.43«+i.sin 2.430. 

Of this the imaginary part, rejecting i, is: 

A-g-*-"f (.0002183 sinL2.43<-. 0001283 cos 2.436). 

Turn next-to the residue, omitting the factor— J, at X= — 4.18— 2.43i, or 

J? .02861 (-4.18-2.43i) f^"«(c03 2.43f-i sin 2.43f) 
(-4.115-2,243i) (-4.115- 2.61 7i)» (-4.86i) 

The angle of the coefficient is 26°.65; the log mag is 6.3964. Hence 

5,= (.OOD2226 + .0001117i) r*-^'* (cos 2.43«-i sm 2.46/). 

Of this the imaginary part, neglecting i, is: 

B=«r«.i«(_. 0002226 sin 2.43t+. 0001117 cgs 2.43iO. 

Adding A and B, the dependence of 6 on the short oscillation is 

-fl/J =<-*■»«(- .0000043 sin 2.43<- .0000168 cos 2.43*). (3) 

The effect is very small, indeed quite negligible (compared with (4) below) as might be imagined 
from the small results f oimd for other types of gust in the two previous papers. 
The residue at X = - .0654 - .187i is 

» .02851(-.0664-.187i) e-'^Cco3 .187f-tsin 187f) 
^ (- .374i)»(4.115 +2.24St) (4.11fe-2.6l7i) 

mag (-.0654-.187t)=».1981, log mag =9.2969, ang = 250°.72. 

Theangleof thecoe£aci6ntis74°.58. The log mag is 7.2377. Hence 

i2,=ff-«««(cos .187«-i sin .187*) (.0004596 + .001666i). 
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Of this the imaginarj part is: 

-.0004596 sin .187* +.001666 cos .187*). 

For —61 J there remains only to calculate the residue at X= — .0654+187t, which is the 
hardest of all, since this factor is squared. We hare to find 

d r .02851Xg» Jl 

^ L(>- + -0654 + .lS7i) (X + 4.18 + 2.43t) (X + 4.18 - 2.43i) J * 

The dnriTatiTo of a product is often calculated most easily by taking the logarithm before 
differentiation. For: 



The derivative of the logarithm of the bracket wiiih respect to X is 

n 1^, 1 1 1 

■^"X"''' X+.0654+.187i X+418+2.43i X+4.18-2.43i* 

As X = - .0654 + .187t we hare 

^" - .0854 + .187t"^*~T374t ~4.1 15 + 2.617i~4.115 - 2.243i 
= - 1.1629 -4.766t+ *+ 2.674i- .1731 + .llOlt- .1832- .0999* 
= - 1.985-1.9921+* , ang=225.*'10 , log mag =.4491. 

Then D must be multiplied by 

^ .02861(-.0654 + .187i)g-'^«(co8 .187*+tsin .187*) 
^ (.374i) (4.115 + 2.617t) (4.116 -2.243i) 

The angle of the coefficient is 15.°42; the log mag is 6.8106. When multiplied by B we have 

[*(.0008233 + . 000171 6i) + (-.0008935 -.001583)1 «-«««(cos .187*+t sin 187*). 
Taking the imaginary part, we have: 
F= [* (.0006233 sm .187*+ .0001716 cos .187*) + (-.0008935 sin .187*-.001583 cos .187*)]. 

On adding 0 and F the eSect on the long oscillation is 

-fl/J=«-«««[*(.000623 sin .187*+ .000172 cos .187*) 

+ ( - .001363 sin .187*+ .000083 cos .187*)]. (4) 

The accuracy, of coiu:se, at this point is not great. Tlie total result for B should give 0=0 
and d6/dt=0 when *'=0, and it does within the estimated remaining accuracy. 

Although the short oscillation is of importance if the values of dw/dt or dBfdt are desired, 
it is (as seen above and in previous papers) of very little use in considering the variable w or 0 
at least in forced motions where the applied force operates relatively slowly (mild or moderate 
as distinguished from sharp gusts) ; it is quite insignificant for the path. The above calcula- 
tions could therefore be abridged somewhat by the device mentioned in my second paper of a 
pfo'tial resolution with partial fractions: 

1 .016X + .089 -.01601X+ .04263 

(X+ .0664 ± .187*) (X+ 4.18 ± 2.43i) "(X+4.18±2.43»)''' (X+.0654±.187i) ' 

The residues to be cdculated for —O/J are those of 

.02851X(.016X+ .089) , .02861X(-. 0160 1X+. 04263) 

(X+4.18±2.43i) (X+ .0654-. 187*) (X+.0664+.187t) (X+.0654-187i)»' 
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In the first there is no residue for X= —.0654— .187* and the residue for X= — .0654 +.1871 is 
negligible compared with that of the second expression for \= — .0664 + .187i (which is the easy 
one to calculate). As the residues at X^^ — 4.18±2>43i refer to the short oscillation, they 
need not-be calculated, and therefore the only calculations really necessary for discussing 
the path are those for the residues of the second expression. The work was carried through 
on this basis and checked with that obtained by the complete analysis above to witbin 2 or 3% 
On the abbreviated plan just outlined let us calculate the result for — w/J. The residues 
are required for : 

X*(.557X+2.458)e^'(-.01601X+.04263) 
(X+ .0654+.187i) (X+ .0654 - .187i)' 

at X= - .0654- .187i and X= - .0654 + .187i. The firet is (- .0654 - .187i)» (2.422 - .1042i) 

(.04368 + . 0029941) = <r-''«« (cos .187*- 1 sin .187*) (.02375-. 01 797i). 
Of this the imaginary part is 



Gt= «-.<««( -.02375 sin .187*- .01797 cos .1870. 

The second is 

^ p(.557X+ 2.458) e'^'(-.01601X+.04263) "| 
^ X+. 0654 + .187* Jx=-.0654+.187t. 



Apply logarithmic differentiation as before, and the residue is seen to be the value of the 
bracket multiplied by 

2 .557 .01601 1 

X"^.567X+ 2.468'^ -.01601X + . 04263 X + .0654 + .187?' 

or 

-.065U.187j + 2.422-!l042* +<- .0436fZ99 -4;-:3ls- -3-258 -9.532i+* + 2.674i + .2295 

+ .009871 - .3648 - .02506i =*- 3.391 - 6.873t 

The value of the bracket itself, apart from , is - .006718 + .008884t. 
Hence the residue is 

<r'«««(cos .187«+i8in .187*) [<(-. 006718 + .0.08884i) + (.08385 + .01607i)]. 

Of tiu3 tlifl imaginary part is 

006718 sin .187*+.0160 cos .187*) + (.08385 sin .187*+. 01607 cos .187*)]. 

Adding G and H, the result for— w/J is 

-w//-r^''«[*(-. 00672 sin .187*+. 00888 cos .187*) + (.0610 sin .187*-. 0019 cos .187*). 

The last term should, of course, check out so as to give zero when <=0, but owing to the 
omission of the terms corresponding to the short oscillation the check can not be expected to 
be exact. Moreover the derivative should also vanish, but has the value + .02. This would 
correspond to a term —w/J =• .005e~^-**cos 2.43*. Now when treating the gust .7(1 - «^0, which 
in its initial effects should not be far different from e/g--'"* sin .187*, the term — w//=e~^ "« 
(.004 cos 2.43*— ,003 sin 2.43*), for which the derivative has the value —.024, was actually 
found. Such a failure to check as occurs in the value of —w/J here determined can not be 
regarded as anJndication of error in the calculation; and an independent calculation checks well 
with the value above given for —w/J. 

The path in space is given by 

a- f(.v>+n5M)dt=J r«r-«««[*(-.062 8in .187*-.029 cos .187*) + (.095 sin .187*- 

.008 cos .187*)](i<. 
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To integrate an expression Uke e~^t cos ht or e~<^i sin ht, the simplest thing is to integrate 
te^ "ie^'**^^, which, may be found in integral tables, and take real and imaginary parts. 

Jie-^^f-^- J-^3^-^3;^^](cos W+isin ht), 

Cs m nfiQ) sia if-"^ cos 6*) , (^-a*) cos U+2ab sin Ml 
jter-^ cos 6i(Zt-«-»|_^ '-+ ^-(^^pp J' 

V, ni ' ,^-*(osin W+& COB if) , (6»-o») sin Jf-2o& COS &n 

J =n_ ra^ w+w ^J* 

Then 

g=/r--»««[f(.34 cos .187i-.03 sin .187t)-1.8 sin .187*]. 

The values for —^J are, for i = 10, 12, 14, respectiTely, about 1.4, 1.9, 2.0; and from then 
on the values decrease. Compare this for —e/J in the case of the periodic gust Ui=J aia .2t, 
where the corresponding values are 1.9, 2.6, 3.0. It is seen that the damped periodic gust is 
decidedly less effective than the undamped gust. This is precisely what I foresaw, and indeed 
what must be admitted a priori unless astonishing powers of discrimination are given to the 
machine. A gust Je"-""* sin .187f or J sin .2t does not differ in. general character from 
the gust of the form J(l— <~-*) or J(l— «~-"'*) during the first rise from zero to a maxi- 
mum — ^the plot of the int-ensitj is nearly a straight line until the maximum is approached and 
the slopes of the lines are nearly equal. The machine is by its inertia an integrating, rather 
than a differentiating device, and should give similar displacements in all four ceses. The 
damping in the first case tends merely to diminish the effect. The maximum rise in the third 
and fourth types is about —S.SJ, in the second case about — 3J (as is natural since the third and 
fourth gusts persist where the second begrus to decrease after about 8 seconds), and His damped 
gust gives — 2J (as is again natiu-al, since the mayimiiTn of that gust is only about .6 J owing to the 
damping factor). The successive forced oscillations in the case of the damped gust drop off 
very rapidly, whereas the straight periodic gust brings a decided resonance into play after the 
natural motion subsides. The conclusion is that the constants of integration are such as to 
ma.sk the effect of the damped gust in the first quarter period, whereas the damping makes the 
effect small at the subsequent times of maxima. 

Defastuent of Petsics, 

MASSACHnsETrs Institdte of Technoloot, 

October IS, 1918. 



